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SUMMARY 

Two methods are presented for determining spectra from filter spectroradiometer 
data. The first method uses only data from the filter spectroradiometer, while the sec- 
ond combines these data with monochromator radiometer data. The result, in either 
case, is a computed spectral curve that is consistent with the experimental data. To il- 
lustrate the data reduction techniques, sample spectra from a carbon arc solar simulator 
are given. The experiment is performed to determine the degree of agreement of mea- 
surements on the same source by the monochromator and filter instruments. 


INTRODUCTION 

Measurement of the spectral irradiance of solar simulator systems is a current field 
of interest. The solar simulator, as a part of a space environment facility, is intended 
to duplicate the solar spectrum. Because the spectrum is only approximated, the differ- 
ence must be found and evaluated to avoid error. 

In optical work, the filter method is of interest because the compactness and sim- 
plicity of the equipment allow field measurements, as opposed to laboratory measure- 
ments, to be made with relative ease. In addition, this method can provide an indepen- 
dent check of measurements made by a prism or grating monochromator. A monochro- 
mator is large and cumbersome, and the fact that the geometry and spectrum of the 
source being measured is different from the calibrated standard lamp used for compari- 
son also causes difficulties. A summary of the difficulties encountered when the mono- 
chromator method is used may be found in reference 1. This situation motivated 
Wyszecki (ref. 2) to use filters for determining the spectral sensitivity of a photoelectric 
detector. The same method also reported by references 3 and 4 is applied to spectro- 
radiometry in this report. 

The filter factor method of data reduction described in the open literature (e. g. , 
ref. 5) was felt to be inadequate. The statements of the method found there are incom- 
plete in that they do not explain how to determine the spectrum of an unknown source. 

The filter factor method requires that the bandwidth limits of a filter be defined arbi- 
trarily. 


The problem presented herein is confined to determining a best estimate of the ac- 
tual spectrum, based on the data available, within the framework of the arbitrary 
choices that have to be made in selecting a reduction method. No consideration is given 
to determining the number of filters and the transmission characteristics required for 
measuring a spectrum. 

Two principal methods of data reduction are presented. Each requires a large 
amount of computation, which necessitates the use of a computer. Neither method gives 
a truly unique solution; hence, any number of equally valid solutions could be found by 
imposing different arbitrary conditions on the solution. 

The first method of data reduction is the '’linear equivalent spectrum” method. It 
is so named because the calculated spectrum, which consists of a series of straight line 
segments, is equivalent to the actual spectrum in the sense that it would give exactly the 
same experimental results for the filter measurements as would the actual spectrum, if 
such an experiment could be performed. The resolution of this spectrum will be limited 
by the number of filters used. This method treats filter spectroradiometry as a com- 
pletely independent method of measurement and utilizes no prior knowledge about the 
spectrum in the data reduction process. It does not, however, rule out the possibility 
that the choice of filters may have been made with some knowledge of the spectrum to be 
measured. 

The second or "composite spectrum” method uses both filter measurements and 
high resolution spectral measurements obtained by means of a monochromator; hence, 
this method utilizes prior knowledge about the spectrum. The procedure is to apply a 
polynomial correction to the monochromator data in order to achieve a best fit of the 
filter data. The final spectrum has the same resolution as the monochromator data. 

Spectra irradiance measurements from a carbon arc solar simulator are presented 
for the purpose of illustrating the data reduction methods introduced. 


ANALYSIS 

The required equipment for filter spectroradiometry consists of the following: (1) a 

set of narrow band-pass filters, each covering a different spectrum region (regions may 

or may not overlap); (2) a detector to receive the radiation passed by the filter; and (3) 

an amplifier and voltmeter to measure the electrical output of the detector. 

th 

The output voltage of the spectroradiometer with the i filter in position is 
given by 



H x (X)T.S(X)dX 


( 1 ) 
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(Symbols are defined in appendix A. ) The R. (A) is the irradiance produced by the source 

An 

at the point of measurement (may be expressed in mW/(cm )(p)) and is a function of the 

•Hi 

wavelength A. The T^ is the transmission of the 1 filter, and S(A) is the sensitivity 
of the detector. If the detector is a blackened thermopile, the sensitivity is assumed to 
be constant with wavelength. Equation (1) also assumes the detector output to be linear 
with flux level. If these are not the case, the experimental voltages must be corrected to 
give values that are proportional to spectral sensitivity and flux level in order to apply 
the theory that follows. 

The data reduction problem is that the Vj are known but that the H^(A) is not. The 
method of this analysis is to construct an assumed irradiance function a (A) such that 
when substituted for H^(A) in equation (1) an output voltage V a ^ results which agrees 
with the V. value obtained experimentally. The R (A) must satisfy all nV. values 

I Aj Q/ X 

simultaneously, where n is the number of filters. This report presents two general 
methods of accomplishing this. The first, an independent method, utilizes only the data 
from the filter spectroradiometer. This method will result ina H. (A) curve that will, 

A, £1 

in general, lack the detailed structure of the actual H(A) curve. The second, a dependent 
method, combines filter and monochromator data. The filters provide a polynomial cor- 
rection term for the monochromator result. The R (A) curve then has the same reso- 

A , el 

lution as the monochromator curve. 

Linear Equivalent Spectrum Method 

The assumed spectrum a (A) consists of n points, one for each filter. Each 
point is defined by the pair of values R „(A) and A.. The wavelength value A. chosen 

A y cX X X 

for each point is the center of gravity of the corresponding filter’s transmission curve. 

It is defined by 



Each irradiance value R, „ . is taken to be unity as a starting assumption, although any 

a, a, 1 

aribtrary set of values, each greater than zero, would serve as well. The spectrum 

H x a (A) is thus defined by the series of straight line segments connecting the n assumed 

points (and extending beyond them in the case of the first and last points). The spectrum 

H-. _(A) thus consists of n - 1 regions. In each region it is defined by a straight line: 
a , a 
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(3) 


a > i+1 a > 1 + jj. 


A, a 


X i+1 ' X i 


A, a, i 


In the event equation (3) would require a negative value, a (X) is taken to be zero. A 
region, of course, extends from the center of gravity of one filter's transmission curve 
to the center of gravity of the next, except for the first and last regions, which extend in- 
definitely beyond the first and last filters. 

The iteration scheme for the solution H. JX) is now given. The starting assumption 

Aj ci 

is 


for each n. This expression defines the trial spectrum a (X) everywhere in accord- 
ance with equation (3). The n values, each a calculated voltage from the assumed irra- 
diance, can be found by solving 


v a ,i=jf' H x,a< x > T l< x > S(X)<iX 

The new value assumed by each H x a i is 

V. 

( H >-, a, i) new y . ( H X,a,i) old 
a, i 

The process is repeated from equation (5) on until 


(5) 


( 6 ) 


- 1 


< 0. 0001 


a, i 


(7) 


for each and every n. A limit of 100 iterations was permitted. In most cases conver- 
gence is obtained within this arbitrary limit. When convergence is completed, the as- 
sumed irradiance function H. 0 (X) is the solution to the problem, and when it is used in 
place of H^(A) in equation (1) it is equivalent in the sense that it gives the same V i val- 
ues as were obtained experimentally. 

Once a measurement is performed and the data reduced, the n ratios can be found: 
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( 8 ) 



H X, a, i 
V i 


On any subsequent measurement of a similar radiation source, an approximate solution 


H A, a, i = R i V i 


(9) 


can quickly be obtained. 

An example of the spectrum produced by this method is the solid line curve of fig- 
ure 1. This result was obtained on a carbon arc solar simulator. The instrument used 
was the Eppley Laboratory’s Mark IV Radiometer equipped with a set of 12 filters, 
which the manufacturer chose for carbon arc spectra and designated as "C” filters. 


Applications of Linear Equivalent Spectrum Method 

It is now possible to compare the linear equivalent spectrum to a higher resolution 
spectrum such as might be obtained on a monochromator. The comparison can be made 
by converting the monochromator result to what the filter method would produce. The 
V i is calculated according to equation (1) where the H^(A) is the actual high resolution 
monochromator spectrum. The calculation proceeds as though the V. had been obtained 



Figure 1. - Spectroradiometric measurements of carbon arc solar simulator. 
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Figure 2. - Grating monochromator measurement of carbon arc solar simulator. 

experimentally with a filter radiometer, and the equivalent spectrum is derived by equa- 
tions (4) to (7). This equivalent filter spectrum would have been obtained by measuring 
an actual source with the spectrum indicated by the monochromator. 

An an example of the conversion of a high resolution spectrum by use of the equiva- 
lent spectrum method form, figure 2 shows a grating monochromator measurement of the 
same carbon arc solar simulator. The measurement was made in the same location as 
the filter measurement. The bandwidth of the instrument was as low as 0. 0008 micron 
in the 0. 31 to 0. 7 micron range. The data were taken at 0. 01-micron intervals. The 
transmission curves of the same set of 12 filters of the Eppley Laboratory's Mark IV 
radiometer were used with this spectrum. 

The equivalent spectrum values at the appropriate center wavelengths are plotted in 
figure 2. Linear interpolation between these values generates the equivalent spectrum. 
Both the high resolution spectrum and its equivalent spectrum yield the same detector 
voltages if calculated by equation (1). Since the same set of filters were used to calcu- 
late the equivalent spectrum of the monochromator measurement and for the actual filter 
measurements, direct comparison of the two instruments is possible. 

The equivalent spectrum of the monochromator is plotted as a dashed line in figure 1 
for comparison to the filter measurement. The difference between these two curves in- 
dicates possible accuracy of the instruments and can give a better estimate of the true 
spectrum. 
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120 . 


Linear equivalent form of 

filter measurement 



Wavelength, [i 

Figure 3. - Comparison of filter measurement with zero air mass solar irradiance. 

A second attractive application of the equivalent spectrum method is to determine the 
difference between the zero air mass spectrum of solar irradiance (Johnson spectrum, 
ref. 6) and the spectrum of a solar simulator. The equivalent spectrum of the Johnson 
spectrum was calculated as previously described for the monochromator data. It is plot- 
ted in figure 3 and compared to the filter radiometer equivalent spectrum of the carbon 
arc simulator shown in figure 1 (p. 5). The difference of the curves is the deviation of 
the simulator from the ideal spectrum. The equivalent of the Johnson curve was plotted 
with its peak at 100. The area of all curves illustrated in this report is equal to the area 
under the Johnson curve between 0. 25 and 2. 2 microns. 


Composite Spectrum Method 

The linear equivalent spectrum method made use of no prior knowledge about the 
spectrum. A composite reduction method has also been developed that utilizes the strong 
points of both monochromator and filter radiometer measurements. (Differences between 
the monochromator and filter radiometer results are shown in fig. 1 . ) The monochroma- 
tor is quite reliable for giving the relative intensities of two spectral regions close to- 
gether in wavelength. The characteristics of the instrument vary, however, as the wave- 
length is changed so that there is less confidence in the relative intensities as the wave- 
length becomes more separated. The filter radiometer because of its simplicity is less 
prone to this type of error. Thus it appears attractive to use the monochromator for its 
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high resolution content and supplement it with filter measurements. A method generally 
applicable is to seek a solution of the form 


k-1 

= a®) 

j=0 

where k < n and m (X) is a high resolution spectral curve resulting from a monochro- 
mator measurement. The set of Pj values is then found that gives the best agreement 
between the calculated V . from equation (5) and the experimental V- indicated in 

M-J J. 1 

equation (1) . This method has the effect of using the monochromator to define the spec- 
trum over the narrow range and the filters to define broad spectrum regions with respect 
to one another. The filter measurements determine the polynomial correction term to be 
applied to the monochromator data. 

It is possible to find the solution of the form indicated in equation (10) by an analytical 
least squares method, which requires no iterative process. This is shown in appendix B. 
As a special case of this procedure, if EL m = 1, a polynomial solution for EL „ can be 

A, m cl 

obtained that makes no use of prior knowledge of the spectrum. 

The monochromator and filter measurements of figures 1 and 2 were used as an ex- 
ample of the composite method. A polynomial of second order (p Q + PjX + P 2 X 2 ) was 
used in equation (10). The value of e (appendix B) was chosen as 1 percent of the maxi- 
mum V. obtained experimentally. The calculation yielded a composite spectrum of high 
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TABLE I. - COMPARISON OF MONOCHROMATOR 


AND FILTER MEASUREMENTS 


Center 

Deviation of 

Deviation of 

wavelength 

composite 

uncorrected 

of filter 

spectrum using 

monochr omator 


second-degree 

spectrum, 


polynomial, 

percent 


percent 
(fig. 4) 

(fig- 1) 

0.264 

-2.8 

-1.7 

.320 

10.9 

12.7 

.355 

-10.7 

-9.0 

.406 

-.3 

1.9 

. 459 

-1.8 

.7 

. 548 

3. 9 

6.7 

.643 

3.8 

6.7 

.753 

-6.3 

-3. 9 

.880 

-6.0 

-4.3 

1. 167 

5.4 

4.6 

1.543 

.7 

-6.0 

2.082 

-1.2 

-16.4 

Root- mean- square 

5.6 

7.6 

deviation 




resolution as in figure 2, but the correc- 
tion applied by the filter data caused a low 
resolution slowly varying deviation. In 
order to compare the results of the com- 
posite method with the filter radiometer 
alone, the composite spectrum was con- 
verted by use of the equivalent spectrum 
method previously applied in figure 2. 

The equivalent spectrum by the com- 
posite method is plotted as a dashed line 
in figure 4 for comparison to the equiva- 
lent spectrum of the filter radiometer 
alone. The deviations, at the center wave- 
length points, between the monochromator 
measured curves of figures 4 and 1 and 
the filter radiometer measured spectrum 
are tabulated in table I. The difference of 
the deviations in the two columns when 
plotted against wavelength form a second- 
degree polynomial. The root-mean- 
square deviation of the composite spec- 
trum has been reduced from 7. 6 to 5. 6 


percent. The correction is only partial with the second-degree polynomial. 

The biggest problem in using this composite method is obtaining a good monochroma- 
tor curve to use with the filter data. The two sets of data should be obtained under as 


nearly identical conditions as possible. It is not proper, for example, to take a mono- 
chromator curve on a bare carbon arc and combine this with filter data taken on a solar 


simulator. The spectrum is modified in passing through the optical system of the simu- 
lator, but even more important is the fact that the simulator may utilize only a portion of 
the source or perhaps a different portion of the source than is seen by the monochroma- 
tor in making its measurement on the bare arc. 


CONCLUSIONS 

The first data reduction method - the linear equivalent spectrum method - yields a 
spectrum that uses only voltage outputs of a detector and the transmission characteris- 
tics of the filters. No prior knowledge of the source is required. The method is mathe- 
matically consistent in the sense that the computed spectrum yields the same radiometer 
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voltage readings as were obtained by measuring the incident radiation. This condition is 
necessary but not sufficient to obtain the best estimate of the actual spectrum with a given 
set of filters. Two other necessary conditions not investigated in this report are being 
investigated at Lewis Research Center. One involves the optimum method of interpola- 
tion. The other is the choice of filter shape and bandwidth which affect the uniqueness of 
the spectrum in the sense that a given set of filters if displaced equally in wavelength will 
produce a different spectrum. 

It was shown how to compare high resolution spectra of a monochromator with the 
low resolution of a filter radiometer. The comparison was applied to determine the devi- 
ation of measurements of the two instruments from each other. It also enabled a better 
estimate of the spectrum of a source, and led to the composite method of superimposing 
the high resolution of a monochromator spectrum on the filter radiometer spectrum. The 
accuracy of the filter radiometer and the resolution of the monochromator were retained 
during this comparison. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 22, 1965. 
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APPENDIX A 
SYMBOLS 


G 

error term 

€ 

number 

H X 

irradiance (spectral) 

X 

wavelength 

V 

s+OO 

/ H. T,SX ] dX 

Subscripts: 

J 0 X, m 1 

a 

assumed 

p i 

coefficient in monochromator 
correction polynomial 

i 

integer referring to particular 
filter 

R i 

ratio of assumed irradiance to 
filter detector voltage 

j 

integer referring to term in mono- 
chromator correction polynomial 

S 

detector sensitivity 

k 

number of terms in monochroma- 

T i 

filter transmission 


tor correction polynomial 

V 

detector output voltage 

m 

monochromator 

x i 

V- - V • 
1 a, l 

n 

number of filters 
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APPENDIX B 


ANALYTICAL LEAST SQUARES FIT OF FILTER DATA 


It is desired to find a solution of the form 


p j xi 

3=0 


where k ^ n. The H must agree with the data from the n filters; that is, the 

SL 


v a,i=/“ H X,a< X > T i« S < X > dX 


must agree with the V.. measured experimentally: 


V a,r /'” H X,a< x > T i( x ) S ( X ) dX '/" H X,m( P 0 + P l X+ ' ' ' +p k-l X ‘ I " 1 ) T i S 


V a,rX” Hx >“ PoTiSdX+ X"° Hx - mPlXTiSdX + ' ' ' + /" Hx ,m P k-l xk ' lT i Sdx 


s+OO /* OO 

V n . = P n / H. T.S dX +P 1 / H. T.S\dX 
a, 1 0 J Q X, m i 1 J 0 i 


4 * • • * 4 * 


P k-1 M B1 > 


where in general 


Then, 


V a,i " P 0 K 0,i + P l K l,i + * ‘ * +P k-l K k-l,i 


i,i' /'” H X,m( X > T i« S < X > xidX 


k-1 

V n f = Y PaK- . 
a, i Lj ] ],i 

3=0 
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with the coefficients Kq Kj j . . 

. K. , . evaluated once by numerical integration. 

Let 

7 

v. - V 

^ i a, 1 

(B2) 


1 V i +e 


and define 

n 




(B3) 


i=l 



If e = 0, X. is the relative error of V ■ compared to V.. Some small value for e is 

1 Clj 1 1 

necessary in the event one of the Vj is zero. The G is the sum of the squares of all 
the individual errors. 

The smaller G is, the better the H-. solution. The problem is to find that set of 

a, a 

Pj values which makes G a minimum. The term G constitutes a merit function that 
tells how good any trial solution is. Once one has taken a measurement, obtained the 
monochromator data m , and arbitrarily selected e and the number of Pj's the solu- 
tion is to have, G is found to be a function of only the P^'s, or G(Pq, Pj . . . P^_ j). A 
necessary condition for G to be a minimum is that 


3G 3G 



dG 


BP- 


k-1 


= 0 


(B4) 


ap o 9p o L^i V v i +e 

i=l i=l 


(B5) 


n 


_0G 

0P 


'r-+Y ^-2 (V i- P 0 K 0,: 

o Z— /(V i+e ) 2 

i=l 


p i K i,i 




Rearranging the expression gives 
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n n n n 

y Xi K o.i P v vm ,,„ p 

} /—J(V i+ e ) 2 1 Z_^/(V i +e) 2 k_1 Z_-/ (V i+ e ) 2 Z_-/(V i+ e) : 

i=l 1 i=l 1 i=l 1 i=l 


Similar expressions can be found for 9G/3Pj, SG/SPg, and so forth. Each is then set 
equal to zero, which gives 



These k simultaneous linear equations in k unknowns can then be solved for the P^'s, 
which completes the solution. 

When this method is applied, judgment must be exercised in choosing the number of 
p.»s used. In no case can it exceed the number of filters. The fewer P^'s used, the 
more the data is smoothed and the less resolution obtained in the correction term. In 
practice, the number of Pj's that may be used is limited by round-off errors in the 
computation process. 

The choice of e determines the type of merit function. If e is much smaller than 
the largest value, equal weights are given to equal percentage errors in matching 
the computed V„ • to the experimental V-. In this case the e merely prevents division 
by zero in the event a Vj is zero. But if e is much larger than the largest V., equal 
weights are given to equal absolute differences between the and V a .. Here again 
judgment must be exercised. 
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